Passage Times for Unbiased Polymer Translocation through a Narrow Pore 
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We study the translocation process of a polymer in the absence of external fields for various pore 
diameters b and membrane thickness L. The polymer performs Rouse and reptation dynamics. The 
mean translocation time (rj) that the polymer needs to escape from a cell, and the mean dwell time 
{rd) that the polymer spends in the pore during the translocation process, obey scaling relations in 
terms of the polymer length N, L and b/R g , where R g is the radius of gyration for the polymer. 
We explain these relations using simple arguments based on polymer dynamics and the equilibrium 
properties of polymers. 

PACS numbers: 36.20.-r, 82.35.Lr, 87.15.Aa 
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Transport of molecules across membranes is an essen- 
tial mechanism for life processes. These molecules are of- 
ten long, and the narrow pores in the membranes do not 
allow them to pass through as a single unit. They have 
to thus squeeze — i.e., translocate — themselves through 
the pores. DNA, RNA and proteins are such naturally 
occuring long molecules Q, 0, 0, in a variety of bi- 
ological processes. Translocation is used in gene therapy 
0,0, and in delivery of drug molecules to their activation 
sites [3, ISl • Understandably, the process of translocation 
has been an active topic of current research: not only 
because it is a cornerstone of many biological processes, 
but also due to its relevance for practical applications. 

Translocation is a complicated process in living organ- 
isms — the presence of chaperon molecules, pH, chem- 
ical potential gradients, and assisting molecular motors 
strongly influence its dynamics. Consequently, the trans- 
location process has been empirically studied in great va- 
riety in biological literature |l(J, . Study of transloca- 
tion as a biophysical process is more recent. Herein, the 
polymer is simplified to a sequentially connected string 
of N monomers as it passes through a narrow pore on 
a membrane. The quantities of interest are the typical 
time scale for the polymer to leave a confining cell (the 
"escape of a polymer from a vesicle" time scale) [l2j . 
and the typical time scale the polymer spends in the pore 
(the "dwell" time scale) [13| as a function of N and other 
parameters like membrane thickness, membrane adsorp- 
tion, electrochemical potential gradient, etc. [T^. 

These quantities have been measured directly in nu- 
merous experiments |l5j| . A number of (mean-field type) 
theories have be en p roposed for the scaling of these typi- 
cal times [12) El EH during the last decade as well. They 
describe translocation as a first-passage or Kramer's 
problem over an entropic barrier in terms of the "reac- 
tion coordinate" m alone. Here m is the number of the 
monomer threaded into the pore (m = 1, ...,N), and 
the transition rates from m to m ± 1 and vice versa are 
obtained from the derivatives of the free energy w.r.t. m. 

In the context of unbiased polymer translocation (i.e., 
in the absence of external driving fields), the prediction 
of the mean-field theories (which only consider polymers 
with simple random walk statistics) is that the dwell 



time scales as N 3 for Rouse dynamics and as A 2 5 for 
Zimm dynamics |13|. These theories indeed provide in- 
sight into the process of translocation, but their usage 
of the (equilibrium) free energy to determine the transi- 
tion rates from m to m ± 1 implicitly assumes that at 
a fixed reaction coordinate m the polymers equilibriate 
much faster than the typical time for the reaction coor- 
dinate to change its value by ±1. This is not necessarily 
the case for longer polymers or polymers in higher spa- 
tial dimensions [ly . Verifying the scaling results of these 
theories using simulations, too, remains a computation- 
ally significant challenge since it involves simulating long 
polymers and correspondingly long time scales. 

The purpose of this letter is (a) to report the results 
of extensive lattice-based Monte Carlo simulations of the 
unbiased translocation process in three spatial dimen- 
sions, for a variety of polymer lengths, pore diameters 
and membrane thickness; and (b) to trace the physical 
origin of their differences from the existing mean-field 
theory results. Our system consists of two cells A and 
B, each of volume V, that are connected by a pore of 
diameter b in a membrane of thickness L. The polymer 
is modeled as a lattice polymer of N monomers, obeying 
self-avoiding walk statistics. Its movement consists of sin- 
gle monomer jumps to neighboring lattice sites. Jumps 
along the contour of the polymer, i.e., reptation moves, 
are attempted with a higher frequency than jumps that 
displace the contour of the polymer laterally to cause 
Rouse dynamics. A detailed description of this lattice 
polymer model, its computationally efficient implemen- 
tation and a study of some of its properties and appli- 
cations can be found in Ref. 0, . Hydrodynamical 
interactions are not incorporated in this model. 

In our simulations, the polymer repeatedly moves back 
and forth from one cell to the other through the pore 
[see Fig. ^a)]. Our primary interest lies in the scaling 
behaviour of two quantities, (i) the mean translocation 
time (rt), the time required for the whole polymer to es- 
cape from one cell to the other, and (ii) the mean dwell 
time (rd), the time that the polymer spends in the pore 
during the translocation process. To define these quan- 
tities precisely, we introduce the following states of the 
polymer. In states A and B, the entire polymer is lo- 
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FIG. 1: (a) Our system with polymers in different states, (i) state A: all monomers are in cell A; (ii) state T (threaded): some 
monomers are in cell A and the rest in cell B; (iii) state B: all monomers are in cell B. (b) A typical translocation process of 
the polymer [using the definition of polymer states defined in Fig. Hial] for our system. 



cated in cell A, resp. B. States M and M' are denned 
as the states in which the middle monomer is located 
halfway between both cells. Finally, states T and T' are 
the complementary to the previous states: the polymer 
is threaded, but the middle monomer is not in the middle 
of the pore. The finer distinction between states M and 
T, resp. M' and T' is that in the first case, the polymer 
is on its way from state A to B or vice versa, while in 
the second case it originates in state A or B and returns 
to the same state. The translocation process in our si- 
mulations can then be characterized by the sequence of 
these states in time (see Fig.^l. In this picture, the dwell 
time is the mean time that the polymer spends in states 
M or T, while the translocation time is the mean time 
starting at the first instant the polymer reaches state A 
after leaving state B, until it reaches state B again. 

Having set both Kuhn length of the polymer and 
the lattice spacing to unity, we conjecture that for thin 
membranes (L ~ 1) (r d ) ~ N 1+2u F(b/ R g ), and ver- 
ify it for 6 = 1 using polymer lengths up to N = 
1200. For narrow pores (b ~ 1), we argue and ver- 
ify that log(r t ) ~ L. We also observe that (r t ) ~ 
VN {b/ R g )-^ 1+2v+ ^- 2 ^/ u F{b/ R g ). Here v = 0.588 is 
the growth exponent for self-avoiding walks, and 7 = 
1.1601 and 71 = 0.68 are exponents related to the entropy 
of a polymer in bulk and near a rigid wall respectively, 
and F(£) is a scaling function; it approaches a constant 
for £ — y 00 and behaves ~ £-°- 38 ± - 08 as £ _ > 0. 

Argument for the scaling of (id): For thin membranes 
(L ~ 1), a scaling relation between (r d ), b and N 
can be obtained by the following hypothesis: (r d ) = 
N a F(b/N@) for some a and f3. We expect that for 
( T d)i b/Rg is a relevant dimensionless parameter that 
determines how easily the polymer can squeeze itself 
through the pore, since the polymer can "feel" the pres- 
ence of the pore only if its radius of gyration R g is 
comparable to the pore diameter b. This implies that 
b/Nl 3 ~ b/Rg => = v, as R g = XN U . Moreover, 
from physical grounds, such a scaling hypothesis means 
that the scaling function F(£) should approach a con- 
stant for £ — > 00. Since for very large pores (i.e., 



b ^> R g ) the polymer no longer feels the pore, (rd) 
should be the time taken by the polymer to diffuse a dis- 
tance R g =>• (r d > ~ R 2 g /(2D N ) ~ R 2 g N ~ N 1+2u , where 
Dn ~ A -1 is the diffusion coefficient of the polymer in a 
dilute polymer solution. The last relation, together with 
the scaling hypothesis implies that a = 1 + 2v, and 

(r d ) ~ N 2 ^F (±) . (1) 

Additionally, since the monomers within the pore move 
along the contour of the polymer, i.e., reptate, (rd) should 
be independent of L, as long as L <C N [l9j . 

Relation between (r d ) and (t u ): During the dwelling 
process, the polymer necessarily has to pass state M at 
least once. Due to the spatial symmetry between cells A 
and B, each time sequence as depicted in Fig.^is equally 
probable under exchange of states A and B. Additionally, 
each time sequence from A to B is as likely as its time- 
reversed counterpart. 

To devise a computationally cheaper method to mea- 
sure (rd) using these symmetries, we introduce an addi- 
tional time-scale (t u ), the mean unthreading time, which 
is the average time that either state A or B is reached 
from state M (not excluding possible reoccurrences of 
state M). Due to time symmetry, the mean time passed 
since the polymer last left state A, until it reached state 
M, is as large as the mean time passed since the polymer 
last left state B, until it reached state M. Consequently, 

fa) = 2(t„>. (2) 

For b > 1, we expect a similar relation between (r d ) and 
(r„) to hold. However, we do not have a suitable argu- 
ment for it, since for larger pores, the properties of the 
polymer as a macromolecule start to play a role |20|. 

Relation between (r t ) and (t^) for b — 1: The fraction 
of time spent in states M and T compared to (rt) equals 
( r d)/(rt). The probability that the polymer is threaded 
exactly halfway is an equilibrium property; hence, the 
sum of the probabilities pm and p' M that the polymer is 
in state M or M' can be obtained from the contribution 
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of these states to the total partition sum Z to t ■ The ra- 
tios /m = Pm/pw and /t = Pm/pt are non-equilibrium 
properties, but as we show below, it is possible to esti- 
mate these ratios accurately from targeted simulations. 
With these quantities, the average translocation time can 
be obtained indirectly from the dwell time, using 



{n) 



(r d ) 



h(l + /m) 



(pm +Pm')/m(1 + h) 



(3) 



We verify Eqs. 11131) and cross-check their consistency 
using direct simulations to measure (r 4 ) and targeted si- 
mulations to measure (t u ), pm + Pm' , Im and /p. We use 
the lattice polymer model of Ref. [13 • 

First, we estimate the entropic penalty paid by the 
polymer in state M or M'. The partition sum Zf,(N) 
of a self-avoiding polymer of length TV, anchored at the 
origin of an infinite lattice, is proportional to Z\,(N) ~ 
H N TV 7-1 , with 7 = 1.1601 and /1 is a lattice-dependent 
non-universal constant. If this polymer is restricted to 
the half-space z > 0, the same expression holds, but with 
an adjusted exponent 71 = 0.68 |2l|. The partition sum 
of a polymer of length TV in state M or M', threaded 
through a narrow pore (b = 1) in a thin membrane 
(L = 1) is then given by the product of the partition 
functions of two separate self-avoiding polymers of length 
TV/2, each having one of their ends anchored at a rigid 
wall, as Z t (L = 1, TV/2, N/2) = [^ N / 2 {N /2)^' 1 ] 2 ~ 
Z 6 (iV)iV" 7+27l " :L . Adding a linear scaling with V to 
Zb(N), this ratio is also the equilibrium probability that 
the polymer is in state M or M'. Hence, we obtain 



Pm 



+ Pm'(N) ~ N-^ 1 - 1 /V. 



(4) 



With increasing membrane thickness L, since the lat- 
tice coordination number in the pore is much smaller 
than in the bulk, we have log .^(L, TV) ~ L, at least 
as long as L <C N (neglecting logarithmic terms). 
Stated differently, the entropic barrier encountered by 
the translocating polymer increases linearly with L. This 
results in an exponential increase for (T t ) [see Fig. |3J: 



log(Tt) ~ L. 



(5) 



Next, we perform dynamical simulations to determine 
/m, /t and (t u ). The simulations start with a polymer of 
length N , threaded halfway in the pore, and the polymer 
originates from cell A. We then wait until the polymer 
unthreads. If it unthreads into cell A, the starting config- 
uration is labeled M', while if it unthreads into cell B, it is 
labeled M. We record the ratio of the number of polymers 
unthreading into cell A vs. cell B for polymer lengths 
N = 20,40,50 and 80, and obtained 1.37, 1.28, 1.33 and 
1.27 respectively for this ratio; i.e., cell A is preferred 
above cell B by a small factor, which does not increase 
with N . This asymmetry can be easily explained by the 
fact that during the translocation process from state A 
to state B, the polymer accumulates folded segments on 



the B side and stretched segments on the A side, which 
makes it more prone to go back to state A than to pro- 
ceed to state B. Thus, with c ~ 1.3 -1 , we conclude that 



/m(A0 



(6) 



In the same simulations, we measured /t/(1 + It), 
i.e., the fraction of time that the polymer is in state 
M before unthreading. From the theoretical ratio of 
1, i, N - i) and Z t (L = 1, N/2, N/2), we 



found (and numerically observed as well) that 
f T (N) ~ JV- 1 . 



(7) 
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FIG. 2: (rd) obtained from (t u ) for L = b = 1 and N up to 



1200. Solid line: (r d ) = [O.SSTV -13 + e.g^" 2 ' 4 ]" 1 . When the 
same data are plotted in the (b/R g )-((Td)/N 1+2 ") coordinates, 
the £-°-38±o.08 scaling for _p(£) is recovered for £ 
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For each TV, we combined the unthreading times into 
a histogram. We then obtained (t u )(N) from a fit of the 
long-times tail of this histogram (see Fig.[5J). We found 
that for short polymers, (t u )(N) ~ N 3 , while for long 
polymers, (t u )(N) ~ iV 2 .4±o.05^ The exp i anation f tnis 

is as follows: a polymer translocating from A to B pulls 
on its segments in A as it accumulates folded segments 
in B. The resulting strains can be released by reptation 
moves (along the contour) that initiate only at the ends, 
or by Rouse moves (perpendicular to the contour) that 
take place anywhere on the polymer. For the first mecha- 
nism, the scaling of (t u ) with TV is that of reptation [i.e., 
(t„)(TV) - TV 3 , which is the same as the mean-field the- 
ory result], but since there are of 0(N) more segments 
on the polymer (where Rouse moves can occur) than the 
two ends (where reptation moves initiate), the second one 
dominates for long polymers, giving rise to the crossover 
seen in Fig. [21 The precise location of this crossover 
depends on the details of the experiment/simulation. 

Having combined Eqs. (IL'I Hi and (|6I7|I . we obtain, for 
6=1, (r t ) scaling as a function of polymer length as 



(r t ) ~ T/jy2+2i/+7-2 7 i+0.22±0.05 



(8) 
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For b > 1 and L = 1, a scaling relation for (r t ) can 
be obtained if we assume that (r t ) is related to (r^) in 
the same way as in Eq. 0. In this case we expect 
the entropic penalty paid by the polymer in state M 
or M' to behave as a function of b/R g — as explained 
above Eq. JQ), this is the quantity that determines to 
what extent the polymer "feels" the presence of the pore, 
i.e., we expect jV _7+27l_1 in Eq. Q to be replaced by 
(N v /6)(-7+27i-i)/^ . Simultaneously, V in Eq. gj) is to 
be replaced by V/b 2 , as the chance for the polymer to 
find the pore increases linearly with the pore area. We 
however expect Eqs. © and Q) to remain unchanged 
as they only concern a threaded polymer. Together with 
Eq. (JJJ, for 6 > 1 and L = 1 this argument leads one to 

( Tt ) ~ VN (b/ 'fl ff )-U+2"+7-27i)/- 'FQi < Rg ) . ( 9 ) 

We performed direct simulations to verify Eq. @ [see 
Fig. • We started with a polymer in cell A, as shown 
in Fig. ^ and recorded the times (for up to 500 different 
runs) it took to reach state B. We then made a histogram 
of these times and deduced (r t ) from its asymptotic slope. 

In summary, we studied unbiased polymer transloca- 
tion for various pore diameters b and membrane thick- 
nesses L, using a lattice polymer model. We found that 
for L = 1, both the mean translocation time (r*) and the 
mean dwell time (tj) obey scaling relations that involve 
functions of b/R g , where R g is the radius of gyration of 
the polymer. We also showed that for b = 1, ln(r f ) ~ L 
and (rd) is independent of L. We explained these results 
using simple arguments based on the polymer's dynami- 
cal and equilibrium properties. Our analysis explains, for 
the first time, how and why deviations from the mean- 



field theory results occur for long polymers. 




FIG. 3: (r t ) vs. b/R a for N = 20, 40, 80; 6 = 1, ... , 31. Inset: 
ln(r t ) ~ L for b = 1 (circles: N = 20, diamonds: N = 40). 

The persistence length of a polymer in translocation 
experiments is equivalent to the Kuhn length used here. 
In experiments electric field effects due to the applied 
bias voltage and hydrodynamical effects [12 are always 
present; these we did not consider here. Nevertheless, as 
our scaling results are based on very general grounds, we 
expect the same scaling forms (involving b/R g for L = 1, 
or ln(r t ) ~ L for b = 1) to hold when the hydrodynamical 
effects arc taken into account, albeit with different expo- 
nents. The effects of hydrodynamics and external fields 
on translocation are two major topics of our onging work. 

We thank Profs. Henk van Beijeren and Erich Eisen- 
ricgler for helpful discussions. 
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